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Abstract — We consider transmission of two independent and 
separately encoded sources over a two-user binary-input Gaus- 
sian multiple-access channel. The channel gains are assumed 
to be unknown at the transmitter and the goal is to design an 
encoder-decoder pair that achieves reliable communication for 
all channel gains where this is theoretically possible. We call 
such a system universal with respect to the channel gains. 

Kudekar et al. recently showed that terminated low-density 
parity-check convolutional codes (a.k.a. spatially-coupled low- 
density parity-check ensembles) have belief-propagation thresh- 
olds that approach their maximum a-posteriori thresholds. This 
was proven for binary erasure channels and shown empirically 
for binary memoryless symmetric channels. It was conjectured 
that the principle of spatial coupling is very general and 
the phenomenon of threshold saturation applies to a very 
broad class of graphical models. In this work, we derive an 
area theorem for the joint decoder and empirically show that 
threshold saturation occurs for this problem. As a result, we 
demonstrate near-universal performance for this problem using 
the proposed spatially-coupled coding system. 

Index Terms — Gaussian MAC, LDPC codes, spatial coupling, 
EXIT functions, density evolution, joint decoding, protograph, 
area theorem. 

I. Introduction 

The phenomenon of threshold saturation was introduced 
by Kudekar et al. [1] to explain the impressive performance 
of convolutional LDPC ensembles [2], [3]. They observed 
that the belief-propagation (BP) threshold of a spatially- 
coupled ensemble is very close to the maximum-a-posteriori 
(MAP) threshold of its underlying ensemble; a similar state- 
ment was formulated independently, as a conjecture in [4]. 
This phenomenon has been termed "threshold saturation via 
spatial coupling". Kudekar et al. prove in [1] that threshold 
saturation occurs for the binary erasure channel (BEC) and 
a particular convolutional LDPC ensemble and conjecture 
that the phenomenon of threshold saturation is very general. 
Empirical evidence of this phenomenon has also been shown 
for additive white Gaussian noise (AWGN) channels in [5] 
and for a class of channels with memory (the dicode erasure 
channel) in [6], using EXlT-like curves. It is known that 
the MAP threshold of regular LDPC codes approaches the 
Shannon limit with increasing left degree, while keeping 
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the rate fixed (though such codes a have a vanishing BP 
threshold) [1]. So, spatial coupling appears to provides us 
with a new paradigm to construct capacity approaching codes 
for BMS channels. 

Spatial-coupling is now being applied to more general 
scenarios. The noisy Slepian-Wolf problem was considered 
in [7] and the authors showed that the phenomenon of thresh- 
old saturation extends to multi-terminal problems. Threshold 
saturation was shown for the binary-adder channel with 
erasures in [8] by considering EXIT-like curves. Spatially- 
coupled codes have also been shown to achieve the entire 
rate-equivocation region for the BEC wiretap channel [9]. 
The effect of coupling has also been observed for K- 
satisfiability, graph coloring and the Curie-Weiss model of 
statistical physics in [10], and for compressive sensing in 
[11] by using spatially-coupled measurement matrices. 

The notion of universality with respect to channel param- 
eters was considered in [12, p. 398] and [13] in the context 
of compound channels. Yedla et al. focused on the notion 
of universality with respect to channel parameters for multi- 
terminal problems in [14] and showed that spatially-coupled 
codes are near-universal for the noisy Slepian-Wolf problem 
in [7]. Preliminary results in [7] show that spatial coupling 
increases the threshold for transmission over the binary-input 
Gaussian multiple-access channel (MAC) as well. 

In this paper, we consider the performance of spatially- 
coupled codes for transmission over the 2-user binary-input 
Gaussian MAC and investigate the phenomenon of threshold 
saturation for this problem. The Gaussian MAC has been 
extensively studied in the literature and is defined by 

y = /iWxW + /i[2]xM + iv. (1) 

The system model is shown in Fig. 1. The channel in- 
puts X[^],X[^] 6 {±1} and the variation in channel gains 
h^^\h^'^^ 6 [0,oo) can be explained either by fading or 
by different power constraints for the two users. The noise 
is a zero-mean Gaussian random variable, with fixed 
variance of 1. The capacity region is defined as the set of all 
achievable rate tuples {R^^\ R^^^), given by the equations 

i?w</(xW;y|x[2]) 

</(x[2];y|xW) (2) 




Fig. 1. The Gaussian MAC 

The corner points of the capacity region are known to 
be achievable by combining successive cancellation at the 
decoder with single-user codes [15]. This method can also 
be leveraged to achieve any point on the dominant face by 
time sharing or rate splitting [16]. The problem of designing 
good LDPC degree distributions was studied in [17] using 
density evolution (DE), where the authors design good LDPC 
codes for a few points in the achievable region (in terms of 
rate). Another approach was shown in [18] for the case when 
both users have the same transmit power, using EXIT charts. 

These optimization procedures exploit knowledge of the 
channel gains to design good codes. However, in practical 
scenarios the channel gains cannot be known non-causally 
at the transmitter (for example, a fading channel). So, it is 
desirable to fix the rate pair for transmission and view the 
capacity region in terms of the achievable channel gains for 
that rate pair In other words, the capacity region is the set of 
all channel gains {h^^\h^^^) that are achievable, i.e., satisfy 
(2). We call this region as the MAC achievable channel- 
parameter region (MAC-ACPR), to illustrate that the capacity 
region is defined in terms of achievable channel parameters. 
The achievable channel parameter regions (ACPRs) were 
introduced in [19] as reliable channel regions, in the context 
of communication over parallel channels. As the channel 
gains are not known at the transmitter, it is desirable to use 
codes which can simultaneously achieve the entire MAC- 
ACPR (in order to minimize the outage probability for fading 
channels). 

Coding schemes which can achieve the entire rate region 
are said to be universal. However, LDPC codes optimized 
for a fixed 2-user binary-input Gaussian MAC need not 
perform well for different channel parameters. Fig. 2 shows 
the belief propagation achievable channel parameter region 
(BP-ACPR), computed via density evolution (DE), of an 
LDPC ensemble optimized for the equal power case [20, p. 
311]. Even though the performance is close to capacity for 
the equal power case, the performance is far from optimal for 
asymmetric channel gains. So we need to consider additional 
constraints when optimizing irregular LDPC codes to achieve 
universality. Spatially-coupled codes have been shown (via 
numerical computation of the BP-ACPR) to provide near- 
universal performance for the noisy Slepian-Wolf problem 
[7]. In this paper, we consider spatially-coupled codes as a 
potential candidate to achieve universality for this problem. 

To simplify notation, we assume that both users employ 
different codes of rate R, from the same ensemble. We use 
spatially-coupled codes for this problem and provide numer- 
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Fig. 2. BP-ACPR for an LDPC code optimized for the equal rate, equal 
power case and the MAC-ACPR for rate 1/2. The degree profiles can be 
found in [20, p. 311]. 

ical evidence for the phenomenon of threshold saturation. 
Analogous to the MAP threshold for point-to-point commu- 
nication, one can define a MAP boundary for multi-terminal 
scenarios. An appropriate GEXIT kernel is also defined to 
construct GEXIT curves for this problem. These GEXIT 
curves satisfy the area theorem by definition and provide an 
outer bound to the MAP boundary. The threshold saturation 
phenomenon is observed to occur towards the MAP boundary 
i.e., the BP-ACPR boundary saturates towards the MAP 
boundary. This provides numerical evidence that supports 
the conjecture that the outer bound on the MAP boundary, 
computed via the area theorem, is tight. The MAP boundary 
for the (4, 8) regular LDPC ensemble is observed to be close 
to the boundary of the MAC-ACPR, thereby showing that 
spatially-coupled codes provide near-universal performance 
for this problem. To the knowledge of the authors, this is 
currently the only practical coding scheme that achieves near- 
universal performance for this problem. 

The paper is organized as follows: In Section II, we 
describe DE for the joint iterative decoder and the GEXIT 
curves used to identify the upper bounds on the MAP 
threshold are discussed in Section II-E. We briefly discuss 
spatial coupling in Section III and density evolution for 
spatially-coupled codes in Section III-B. The results are 
presented in Section IV. 

II. Density Evolution and GEXIT Curves 

In this section, we briefly cover some notation and back- 
ground for LDPC codes, the 2-user binary-input Gaussian 
MAC and density evolution. We will then discuss GEXIT 
curves for the joint iterative decoder and use them to compute 
an outer bound on the MAP boundary. 

A. Background 

An LDPC ensemble can be characterized by its degree 
profiles. Based on standard notation [20], we let A(x) = 



Y,i Ajx'"^ be the degree distribution (from an edge per- 
spective) corresponding to the variable nodes and p{x) = 
Y,i Pix"^ "^ be the degree distribution (from an edge per- 
spective) of the parity-check nodes in the decoding graph. 
The coefficient (resp. pi) gives the fraction of edges 
that connect to variable nodes (resp. parity-check nodes) of 
degree i. Likewise, Li is the fraction of variable nodes with 
degree i. The design rate of an LDPC ensemble is given by 



:(A,p) = l 



fa p{x)Ax 
X{x)dx 



B. The joint iterative decoder 

Let X[^],X[^] 6 {±1}" denote the input of users 1 
and 2 respectively and Y denote the received vector We 
consider the case when both users employ LDPC codes for 
transmission. To simplify notation, we assume that both the 
users encode the data using LDPC codes from the standard 
irregular ensemble LDPC(n,A,p) and that the transmission 
is bit-aligned. The factor graph of the joint decoder (see 
Fig. 3) consists of two single user Tanner graphs, whose 
variable nodes are connected through a function node [20, p. 
308]. The variable nodes that are connected via the function 
node are chosen at random'. We also assume that the joint 
decoder iterates by performing one round of decoding for 
user one, followed by one round of decoding for user two. 
Let X, = and X = {X^^\X^^^). Without 

loss of generality, we can label the elements of {±1}^ by 
integers X = {0,1,2,3} using the map tt : X ^ {±1}^^ 
defined by ( + 1,+1),1 ^ ( + 1,-1), 2 (-1,+1) and 
3 1-^ (-1,-1). Let 7ri,7r2 : X {±1} be the projections 
onto the first and second coordinate respectively. Then, the 
canonical representation of the channel output is given by 

•^XiiVi) = PY\Xl^]^X^^]{y^Wl{^l),^^2{x^)) 



: exp 



2ct2 



Let m':-'^ f and to?"'! denote the "variable node to function 
node" and "function node to variable node" messages^, 
respectively, for variable node i of the jlh user Here j e 
{1,2} and i e {l,2,...,n}. The message passing rules at 
the function node are given by 



[2] 



log 



log 



[2] 



(3) 



(4) 



Note that this function node operation is not symmetric with 
respect to the users in general. The operations are the same 
only for the case of symmetric fading coefficients i.e., when 



C. Density evolution 

Density evolution is employed to analyze the asymptotic 
performance of the ensemble LDPC(n,A,p). The lack of 
symmetry means that one cannot use the all-zero codeword 
assumption for this problem. Instead, one may assume that 
both users transmit codewords of type one-half, which occurs 
with high probability (a more thorough discussion can be 
found in [20, p. 296]). In this case, px{x) = 1/4, Vx e X. 
We use the notation asAWGNMA = aBAwGNMA(/i[ii,;i[2]) to 
denote the density of the received random variable Y. Let 
/i^2(-,aBAWGNMA) (respectively /2^i(-, aBAWGNMA)) be the 
density transformation operator corresponding to a message 
from user 1 to user 2 (user 2 to user 1) via the function node. 
More precisely, 

/l^2(a,aBAWGNMA) = X! PC^^) /l2 (a (tTi (x) u) , J/j; (m) ) 
xeX 

/2^i(b,aBAWGNMA) = p{x) f2i{h{Tr2{x)u) , ^^{u)) , 

xeX 

where /i2(-,-) and /2i(-,-) are density transformation oper- 
ators corresponding to (3) and (4). Here, a(u) (respectively 
b{u)) is the density of the messages 
These operators can be computed numerically for discretized 
densities following the procedure outlined in [21]. Let ag 
and be denote the L-density^ of the messages emanating 
from the variable nodes to the check nodes at iteration £, 
corresponding to codes 1 and 2. Then, the DE equations for 
the joint decoder are given by 

a^+i = ,f2^i(^L {p{be)) ,aBAWGNMA) ® Hpi^e)) 

be+i = /i^2(-^(/5(a^)) ,aBAWGNMA) ® A(p(b^)), 

where A(a) = ZiA,a®('-i), L(a) = J:^L^a®^'-^\ and 
p(a) = EiPia'^^'"^^. Here, L{p{-)) is the density of the 
messages from the variable node to the function node. These 
equations accurately represent the evolution of densities at 
the decoder due to the symmetry of the variable and check 
node operations. The fixed points of density evolution are 
the triples (aBAWGNMA, a, b) which satisfy 

a = /2^i(i(p(b)) , aBAWGNMA) ® A(p(a)) 

b=/i^2(i(p(a)) , aBAWGNMA )®A(p(b)). (5) 

Let A+oo denote the delta function at +oo and let h^^^ = 
for some A e [0, +oo]. The BP threshold is defined 

by 

hy^^'^^ {\, p, A) = supjft.'^^^ : The fixed point equation 

(5) has a solution (a,b) + (A+oo, A+o^)}, 

/i[^]^«P(a,p,a) = a;iW'«P(a,p,a). 

The set of all points {h^^^' , Ah^^^' ) 3 h^^^' >h^^^'^^ {X, p, A) 
is called the BP-ACPR and its boundary is called the DE 
boundary. 



Other matching rules result in a different performance in general. 
Here, the messages are in the log-lrkelihood domain. 



^^Since the transmission alphabet is {±1}, the densities are conditioned 
assuming the transmission of a +1. 
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Fig. 3. Tanner graph of the joint decoder. The variable nodes of each code 
are connected thi'ough function nodes, which receives the channel outputs. 
The joint decoder iterates by passing messages between the component 
decoders. 

D. GEXIT curves 

The following approach can be applied to any binary-input 
MAC characterized by a single parameter, whose density is 
differentiable and degraded with respect to that parameter 
A through discussion of channel degradation can be found 
in [20, p. 204]. The binary-input Gaussian MAC defined 
by (1) satisfies these conditions for the parameter a, with 
hi = a and h2 = Aa, for any fixed A e [0, +oo). We note 
that any differentiable parameterization would suffice for the 
following discussion. 

Now, suppose that the ith bit is transmitted through a 
channel with parameter q, and that each a, is a differentiable 
function of some parameter a. We omit the dependence 
of the channel on A throughout this section for clarity of 
notation. The GEXIT curve is defined by 



1 " (9 Bn/ 
n ~{ oai aa 



A more convenient expression for the GEXIT curve can be 
derived following the procedure given in [22] for non-binary 
codes. Let y^, = y\yi, (t>i{y^i) = {px,\-ir.i{x\y~z),x € X} 
and $i = (j)i(Y^i) be the corresponding random variable. 
Note that (j>i{-) is the extrinsic MAP estimator of Xi^. 

Lemma 1: Suppose that all bits are transmitted through 
channel with parameter a. Then, the ith GEXIT function is 
given by 

gi(") = Z / aa;.i(u)Ka;(u)du, 

where a^^ i(u) is the distribution of $i given Xi = x and the 

"^To see this, write 

PY„i|Xi(y~iF'i) = — ^ — PY„i(y~i)= „ PY„,(y~i). 



GEXIT kernel is given by 

d 



. , Zt' u\x']p(y\x') , 

p(y|x)log, \ f ^ ' dy, 

oa uYx\p(%)\x) 



(6) 



where denotes the jth component of u. 

Proof: Suppose that each bit is transmitted through a 
channel with parameter and consider the term 

H{-K\Y) = H{X,\Y) + H{-K.,\X,,Y) 

Note that the second term of the decomposition does not 
depend on the channel at position i. So, we get 

d 



da. 



H{X,\Y„<^,). 



We have. 



X X,p{x ,y,(p) 



Y,p{x) Jp{<l>\x) Jp{y\x)\og2 - 



Y,p{x'\(f>)p{y\x') 



p{x\(j))p{y\x) 



-dy 
/ 



which follows from the fact that 

Px,.,Yi,<i>, {x, y, 0) = p{y\x)p{(t)\x)p{x) , 

since Yi ^ Xi ^ Taking the derivative and noting that 

p{xi\(j>i) =p{xi\y^.i), we obtain' 

d 



giio^t) = T.P(x) J p{(j)\x)^J 



-p{y\x) 



Y,p{x'\(t>)p{y\x') 
p{x\(l))p{y\x) '^yj'^^ 
= Y,P(^) / a:r(u)K2;(u)du, 

u 

and the result follows by setting ai = a. ■ 
The GEXIT function is hard to compute and hence we 
use the BP-GEXIT function instead [23]. The BP-GEXIT 
function is obtained by replacing the MAP extrinsic estimator 
with the corresponding BP estimator. Let ^^^'^-^ denote the 
BP extrinsic estimate of Xi after £ iterations of the joint 
decoder. The BP extrinsic estimate is computed using the 
computation graph of depth £ for function node i. Define 
the BP-GEXIT function at the M iteration g^^'^''^{a) in 
a similar manner to [23] (taking an expectation over all 
possible computation graphs) and the asymptotic BP-GEXIT 
function g^^(a) = lim£^oo limn^oo g^^'^'"(a)- For fixed £, 
in the limit of n ^ oo, the computation graph becomes 
tree-like because the computation graphs of the two variable 
nodes (which themselves become tree-like) connected to the 
function node do not overlap with high probability. The 
extrinsic estimate of Xi can then be computed via the 



where e[xi] is the standard basis vector with a 1 in the Xi-th coordinate 
and use the result in [20, p. 29]. 



''The terms obtained by differentiating with respect to the channel inside 
the log vanish. 



extrinsic estimates of Xp-' and X^-'. The asymptotic BP- 
GEXIT function can be computed through the fixed points 
of density evolution (aBAWGNMA(Q)i b) which satisfy (5) 
and is discussed in the following Lemma. 

Lemma 2: Consider transmission over the channel 
aBAWGNMA(a) and let (aBAwGNMA(a) , a, b) be a fixed point 
of DE. Define the BP-GEXIT value of the fixed point by 

G'*''(aBAWGNMA, a, b) = ^ p{x) / Fj;[a, b](u, v)kx{u, v)Aud.v. 

The GEXIT kernel is defined as in (6) and the 

operator Fa;[-,-] computes the density of the extrinsic BP 
estimate given X = x. The BP-GEXIT curve ^^{a) is 
given in parametric form by (a, G(aBAWGNMA(a) , a, b)). 
Proof: Let = u. Then, 

u = 4>{y~i) = {p{xi\y^i),x^ € X} 
= {p{^^l{x,)\y^,)^p{^:2{x^)\y~^),x^ 6 X}. 



If we define 



p(xW=+l|y.O 



u = log 



[1] 



-i|y~0' 



log 



p{X. 



[2] 



-i|y~0' 



then 



_/ e" e^ e" 1 1 e^ 1 1 \ 
"~ \1 + e" 1 + e" ' 1 + e" 1 + e" ' 1 + e" 1 + e" ' 1 + e" 1 + €") 

= f{u,v). 



(7) 



Let z{u) denote the density of U conditioned on xY'^ = +1 
and b(u) be the density of V conditioned on X| = +1. 
Then, a(-u) is the density of U conditioned on = -1 
and b(-t;) is the density of V conditioned on X^ = 
-1. In the limit 71 ^ oo and taking expectation these 
densities ai'e given by the fixed point (aBAWGNMA(a), a, b). 
Let Fj;[a, b](u, u) be the density of <I>f^ conditioned on 
(xW = ^i(x),Xp] = ^2(x)). Then, 

F^[a, b\{u,v) = a {■ni{x)u)h{-n2{x)v) . 

For example Fo[a,b](M,u) = a(u)b(u), Fi[a, b](ii, w) = 
a(u)b(-t;) and so on. The result follows by the definition 
of the GEXIT curve. The kernels Kx(u, v) are defined in the 
sense of (7). ■ 
It can be shown that the BP-GEXIT function is an upper 
bound on the GEXIT function (see the discussion in [20, p. 
206]). The BP-GEXIT curve for the LDPC(3,6) ensemble 
is shown in Fig. 4, for A=l. 

E. The MAP Boundary 

Let h}f^{a) = a and }^f^{a) = Aa, for all i = 
l,...,n. Consider transmission using codes from the en- 
semble LDPC(7i, A, p). For a fixed A, we define the MAP 
threshold as 



MAP 



(A) = M\a : liminf -E[H{X\Y{a, A))] > 
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Fig. 4. BP-GEXIT curve and an upper bound on the MAP threshold 
(computed using the area theorem), for A = 1, of the (3, 6) regular 
LDPC ensemble. GEXIT curves in literature are typically parameterized 
by the channel entropy and the channels get worse as the entropy increases. 
However, the channel gains are a natural parameterization for this problem 
and the channel gets better by increasing the channel gains. So the GEXIT 
values are negative for this parameterization. 

where the expectation is taken over all codes in the ensemble. 
The set of all points {a', A) 3 a' > a^^{A) form the MAP- 
ACPR and its boundary is called the MAP boundary. By 
definition of the GEXIT function, this gives 



f g{a)da= - [ 

j+oo n J+ 



" dH(X|Y(a)) 

00 da 



da 



n 
2k 



H{X\Y{0)) 



The above equation gives us a procedure to compute the 
MAP threshold, using the GEXIT curve. Let a denote the 
largest positive number such that 



J a 



where r(A, p) is the design rate of the ensemble LDPC(A, p). 
Then the MAP threshold a^^^ < a. The upper bound on the 
MAP threshold is shown in Fig. 4 for the (3,6) LDPC en- 
semble, for A= 1. Using this procedure, we can compute an 
outer bound to the MAP boundary by considering different 
values of A. 

III. Spatial Coupling 

Spatial coupling is best described by the {I, r, L) ensemble 
through a protograph [1], [24]. The protograph structure at 
the joint decoder is shown in Fig. 5 for a LDPC(3,6) base 
code. The protograph is generated as follows: Consider the 
protograph of a (3, 6) regular LDPC code. It has two variable 
nodes of degree 3 and one check node of degree 6. Connect 
both the variable nodes to the variable nodes of another 
protograph via function nodes. The resulting protograph 
represents the joint decoder when both users are using (3,6) 



regular LDPC codes for transmission over the 2-user binary- 
input Gaussian MAC. Place 2L + 1 protographs at positions 
-L, L. Each of the 3 edges of a variable node at position i 
is connected to exactly one check node at position i-l, i, 
for each user. 




2L+ 1 

Fig. 5. Protograph of the joint decoder. Shown above are 2L + 1 copies 
of the protograph of the joint decoder for a (3, 6) regular LDPC code. 
The bottom graph shows the protograph of the joint decoder for the 
corresponding spatially coupled code. 



Although this ensemble is very instructive in understand- 
ing the universality of spatially-coupled codes, the EBP 
curves for this ensemble exhibit wiggles around the MAP 
threshold (similar to the single user channels as discussed 
in [1]). The magnitude of these wiggles appears to remain 
constant with increasing L and their presence implies that 
the BP threshold is smaller than the MAP threshold of the 
underlying ensemble. Therefore, the (3,6,i) ensemble does 
not exhibit the threshold saturation phenomenon exactly. To 
overcome this, we use the {l,r,L,w) ensemble introduced 
in [1] for the remainder of this work. 

A. The {l,r,L,w) ensemble 

The {l,r,L,w) spatially-coupled ensemble can be de- 
scribed as follows: Place AI variable nodes at each posi- 
tion in [-L,L]. The check nodes are placed at positions 
[-L, L + w-l], with check nodes at each position. The 
connections are made as described in [1]. This procedure 
generates a Tanner graph for the {l,r,L,w) ensemble. The 
design rate for the {l,r,L,w) ensemble is shown in [1] to 
be 



:{l,r,L,w) 



2L+1 



Two such graphs (generated by the above procedure) are 
taken and the variable nodes (of each graph) at each position 
are connected by a random (uniform) permutation of size 
M via channel nodes. This procedure ensures that all the 
variable node positions are symmetric and enables us to 



write down the density evolution (DE) equations in a simple 
manner, as described in the following section. 

B. Density evolution of the {l,r,L,w) ensemble 

Let ' and ' denote the average density emitted by 
the variable node at position i, at iteration I, for codes 1 and 
2 respectively. Set af^ = bf^ = for i i [-L,L]. The 

channel densities for codes 1 and 2 are denoted by aBMSC and 
bfiMSC respectively. All the above densities are L-densities 
conditioned on the transmission of the all-zero codeword (see 
Section 11). We consider the parallel schedule for each user 
(as described in [1]) and update the correlation nodes before 
proceeding to the next iteration. Let us define 

H(r-1)\ 

I I ~ ^ I I " ~ 1 



J=0 



|H(r-l) j« 



(1 I w-1 \ 

""^ j=0 V''^ fe=0 / 



, H(r-l) 



The DE equations for the joint spatially-coupled system can 
be written as 

= /2^i (r(b|!t+i' b-fL-i)' ^bawgnma) ® 
o(a^^) - a^^^ ) 

,(e+l) f /t^/ {I) (i) \ \ ^ 

- Jl^2 (aj_^^;^,---,a^^^_J,aBAWGNMAj ® 

for i € [-L, L]. For a further discussion of the DE equations 
for the {l,r,L,w) spatially-coupled ensembles on BMS 
channels, see [5]. Using the notation a = (a_i, a^), the 
fixed points of DE are given by (aBAWGNMA,a, b), which 
satisfy 

a. = /(r(b 

l), aBAWGNMA) ® 

bj = / {T{ai-w+i,---, ai+^,_i), aBAWGNMA) ® 

9{^i-w+l,'",^i+w-l)- (8) 

One can use the procedure outlined in [23], known as fixed- 
entropy DE, to compute both stable and unstable fixed points 
that satisfy (8). 

C. GEXIT curves for the (Z, r, L, w) ensemble 

Define the GEXIT value of a fixed point (aBAWGNMA, a, b) 

by 

1 ^ 

G(aBAWGNMA,a, b) = — S G ( aBAWGNMA, ai, bi). 

1 i=-L 

The BP-GEXIT curve g(a), for a fixed A, is the set of 
points (a, G(aBAWGNMA(Q),a, b)). The resulting curves for 
the spatially-coupled (3,6,16,2) and (3,6,32,4) ensembles 
are shown in Fig. 6 for symmetric channel conditions. 
These curves are very similar to the single user case and 
demonstrate the phenomenon of threshold saturation at the 
joint decoder, for symmetric channel conditions. For channel 
parameters not on the symmetric line, these plots imply 
threshold saturation towards the MAP boundary. 
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Fig. 6. BP-GEXIT curves of the (3, 6, L, w) spatially-coupled LDPC and 
(3, 6) regular LDPC ensembles for transmission over a 2-user binary-input 
Gaussian MAC, with A = 1. Also shown is the upper bound on the MAP 
threshold for the (3,6) regular LDPC ensemble, computed using the area 
theorem. 

IV. Results and Concluding Remarks 

It was shown in [1], that for transmission over a BEC, the 
BP threshold of spatially-coupled ensembles is essentially 
equal to the MAP threshold of the underlying ensemble. This 
was observed numerically for general BMS channels in [5], 
[3]. It was observed numerically for multi-user scenarios in 
[7], [8]. The notion of universality with respect to channel 
parameters is important for multi-terminal problems and has 
been discussed in [7], [14], [25]. 

In this paper, we study the 2-user binary-input Gaussian 
MAC and observe that spatial coupling boosts the BP-ACPR 
of the joint decoder to the MAP-ACPR of the underlying 
ensemble. The BP-ACPR for the scenarios considered in 
this paper are shown in Fig. 7 and these results confirm 
the preliminary results reported in [7]. This figure shows 
that spatially-coupled ensembles are near universal for this 
problem. Based on the observation that regular LDPC codes 
with large left degrees behave like random codes and the 
fact that random codes are universal under MAP decoding, 
we also conjecture that increasing the left degree (keeping 
the rate constant) will push the MAP boundary towards the 
boundary of the MAC-ACPR. An analytic proof of threshold 
saturation remains an open problem. Such a proof would 
essentially show that it is possible to achieve universality 
for the 2-user binary-input Gaussian MAC under iterative 
decoding. 

This work can be extended in a variety of ways. For 
example, it is straightforward to dispense with AWGN and 
compute the ACPRs of any suitably parameterized 2-user 
binary-input MAC. One can also generalize these results 
to TO-user MACs, larger input alphabets, and multiple-input 
multiple-output (MIMO) systems. In these cases, the increase 
in computational complexity makes discretized DE infeasible 
and Monte Carlo methods must be used to evaluate the 
DE and GEXIT functions. We conjecture that threshold 
saturation will continue to occur for all these extensions 
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Fig. 7. BP-ACPR of the (3,6,64,5) and (4,8,64,5) spatially-coupled 
LDPC ensembles for the 2-user binary-input Gaussian MAC. Also shown 
are the BP-ACPRs for the (3,6) and (4,8) regular LDPC ensembles. The 
BP-ACPR of the (4,8,64,5) spatially-coupled LDPC ensemble is very 
close to the MAC-ACPR, demonstrating the near-universal performance of 
spatially-coupled codes. 

and that spatially-coupled codes will achieve near-universal 
performance. 
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